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Abstract 

After some introductory discussion of the definition of Finsler spacetimes and their symmetries, we consider 
a class of spherically symmetric and static Finsler spacetimes which are small perturbations of the Schwarzschild 
spacetime. The deviations from the Schwarzschild spacetime are encoded in three perturbation functions <f>o(r), 
0i (r) and 4>2(r) which have the following interpretations: (j>o perturbs the time function, <pi perturbs the radial 
length measurement and <f>2 introduces a spatial anisotropy which is a genuine Finsler feature. We work out the 
equations of motion for freely falling particles and for light rays, i.e. the timelike and lightlike geodesies, in this 
class of spacetimes, and we discuss the bounds placed on the perturbation functions by observations in the Solar 
system. 

1 Introduction 

Since its discovery almost hundred years ago, general relativity has proven to give a very succesful description of 
our universe. Nonetheless, there are good reasons for investigating gravitational theories that are more general than 
general relativity. There are many theoretical predictions, in particular from quantum gravity ideas, that general 
relativity should be replaced by a more general theory at some scale. In order to confront such theoretical predictions 
with experiments, it is necessary to theoretically study all observable effects of the more general theory. This will tell 
by what sort of future experiments deviations from general relativity could be observed, and to what accuracy general 
relativity is verified by present day observation. The PPN formalism provides a mathematical framework for doing 
so; however, it is restricted to metrical theories in the strict sense, i.e., to theories where the gravitational field is 
described by a pseudo-Riemannian metric tensor of Lorentzian signature, as in general relativity. For other theories, 
no such universal framework exists. 

In this paper we want to investigate Finsler gravity theories, i.e., theories where the pseudo-Riemannian metric 
of general relativity is replaced with a Finsler metric. Finsler metrics are characterised by a Lagrangian function 
that is still homogeneous with respect to the velocities, but not necessarily given by a quadratic form. The most 
important feature that distinguishes a Finsler metric from a pseudo-Riemannian metric is in the fact that it breaks 
spatial isotropy even in "infinitesimally small regions", i.e., mathematically speaking, on the tangent space. It is 
true that up to now there is no observational indication for such an anisotropy. (Note, however, Bogovslovsky's [5] 
attempt to explain apparent violations of the GZK limit of cosmic rays as an effect of a spatial anisotropy.) At the 
position of the Earth, deviations from isotropy are strongly restricted by experiments of the Michelson-Morley type 
P3]. However, this result is based on the assumption that the armlength of the Michelson-Morley interferometer is to 
be determined not with the Finsler metric but with an independent Lorentzian background metric. If one assumes, 
by contrast, that the metric which determines the length of solid bodies shows the same sort of anisotropy as the 
metric that determines the light cones, then a Michelson-Morley-type experiment would give a null result. 

Finsler manifolds have been considered as possible spacetime models by a large number of authors. A fairly 
complete list of the pre-1985 literature can be found in Asanov's book [2|. There are several quite different motivations 
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for considering Finsler manifolds as possible spacetime models. Apart from the aesthetic appeal Finsler geometry has 
for many authors, Finsler spacetimes have been recently suggested as a possible explanation for dark matter [B], and 
they have been used in an attempt for explaining the Pioneer anomaly |15j . (As it has now become clear that the latter 
can be explained as a thermal recoil effect [50], also cf. [5] and [35], this motivation should be considered as obsolete.) 
At a more fundamental level, it has been shown that Finsler geometry naturally comes up in models motivated by 
quantum gravity ideas |10] , in particular in Very Special Relativity [5] and in other theories with violation of Lorentz 
invariance [T5] . 

We take this as the motivation for investigating, in this article, the observational bounds on a spherically symmetric 
and static Finsler perturbation of the standard general relativity model of our Solar system. To that end we consider 
the effect such a Finsler perturbation would have on the motion of freely falling particles and light rays. In contrast 
to the above-mentioned Michelson-Morley-type experiments, we will not need any assumption on the behaviour of 
("rigid") extended bodies under the influence of a Finsler perturbation. 

As the class of all spherically symmetric and static Finsler spacetimes is unmanageable (see Section [5] below) , we 
need a special ansatz. As we want to discuss the motion of particles and of light rays, we need a Finsler spacetime in 
which both timelike and lightlike geodesies are well defined. This is an important issue, because in the literature one 
can find many Finsler spacetimes in which the notion of lightlike geodesies is not well defined. As this fact is glossed 
over in many articles, we discuss it in Section [2J below in some detail. Roughly speaking, three different definitions of 
Finsler spacetimes can be found in the literature: The one most frequently, though often implicitly, used in physics 
texts can be found in Asanov's book [2J; an alternative one is due to Beem [3] and a quite recent one, which is a 
generalisation of Beem's, is due to Pfeifer and Wohlfarth [T7]. As we will outline in Section [2] below, none of them is 
appropriate for our purpose: Asanov's definition does not allow to define lightlike geodesies, while Beem's definition is 
slightly too restrictive to define staticity in the most convenient way; the latter observation is unaffected by Pfeifer and 
Wohlfarth's generalisation. Therefore, we will introduce our own definition of Finsler spacetimes in Section [2] below 
which is a slight generalisation of Beem's definition. On the basis of this definition, we will then consider in Section [3] 
a special class of Finsler spacetimes that are perturbations of the Schwarzschild metric. The perturbations preserve 
spherical symmetry and staticity. In this way we arrive at a formalism that allows us to quantitatively study, in 
Section [5] to [8] hypothetical Finsler deviations from general relativity in the Solar system, not only in terms of effects 
on particles but also on light rays; as in general relativity, our light rays are defined as geodesies whose initial vectors 
lie on a unique light cone that determines the causal structure of spacetime. We believe that such a formalism did 
not exist before. It is true that Roxburgh [2T] set up a PPN formalism for Finsler gravity, cf. Roxburgh and Tavakol 
[22] for related material. This, however, was restricted to the very special case of a Finsler metric whose light cones 
coincide with the light cones of a pseudo-Riemannian metric; thereby any Finsler effect on the lightlike geodesies 
was excluded. There is also work by Aringazin and Asanov [T] [3] on Finsler generalisations of the Schwarzschild 
metric and possible observable effects. However, this is based on Asanov's definition for which lightlike geodesies are 
not defined. Earlier work by Coley [7] is also (implicitly) based on Asanov's definition. More recently, Pfeifer and 
Wohlfarth [18] have considered a certain Finsler perturbation of the linearised Schwarzschild metric; however, their 
ansatz is quite different from ours insofar as it introduces birefringence. 

In analogy to the PPN formalism, our analysis will be purely kincmatical, not using any field equation. Several 
attempts of establishing a Finsler generalisation of Einstein's (vacuum) field equation have been brought forward, see 
Rund and Beare [21] (cf. Asanov [2], pp 110), Rutz [25], and Pfeifer and Wohlfarth [IB]. However, it seems fair to 
say that no generally accepted Finsler version of a field equation exists so far. 

As an aside, we mention that Finsler spacetimes in the sense considered in this paper provide a counter-example 
to the Schiff conjecture. In its original version, brought forward by L. Scruff in 1960 [27], this conjecture said that 
a theory must satisfy Einstein's equivalence principle if it satisfies the weak equivalence principle. In our Finsler 
spacetimes there is a unique timelike geodesic for every timelike initial condition, so the weak equivalence principle 
is satisfied. However, as the theory is not based on a pseudo-Riemannian metric, Einstein's equivalence principle is 
violated. 
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2 Definition of Finsler spacetimes and their symmetries 



Historically, Finsler geometry was first established for positive definite metrics. In this case, which is covered in 
standard text-books such as Rund [23] . a Finsler structure is defined in terms of a function F(x, x) that is positive 
and sufficiently smooth on the set of all tangent vectors (x, x) with x =/= 0, and positively homogeneous of degree one, 
i.e., F(x, kx) — kF(x, x) for all k > 0. The Finsler metric is then introduced as the Hessian 

d 2 £{x,x) 

where £(x,x) — F(x,x) 2 , and it is required that this be positive definite for all x ^ 0. The affinely parametrised 
Finsler geodesies are the solutions to the Euler-Lagrange equations of the Lagrangian £(x, x), 

d f d£\ d£ 

dt\dx^J ~ ~dx^ ' ^ ' 

For applications to spacetime theory one would like to have a Finsler metric of Lorentzian signature, and one would 
like to have timelike, lightlike and spacelike geodesies. This requires a modified definition of Finsler structures where 
the Lagrangian £ is no longer positive (hence not the square of a real-valued function F) on the set of all non-zero 
tangent vectors. Therefore Asanov [2, defines a Finsler structure in terms of a function F(x, x) that has the same 
properties as in the positive definite formalism, but is given only on some subset of the tangent bundle which he 
calls the "admissible vectors". The admissible vectors are to be interpreted as timelike. Typically, such a Finsler 
function F involves the square root of an expression that becomes negative on part of the tangent bundle; vectors 
where this happens are not admissible. In Asanov's formalism, timelike geodesies are well-defined as the solutions 
to the Euler-Lagrange equations of the Lagrangian £ = F 2 (or £ = —F 2 , depending on the choice of signature) 
with admissible initial conditions. However, lightlike geodesies are not well-defined. What one would like to define 
as lightlike vectors are the ones on the boundary of the set of admissible vectors; there, however, the Euler-Lagrange 
equations break down because of zeros in the denominator. This formalism of Asanov, in which the Finsler structure 
is well-behaved only on the timelike vectors, is used in many physics papers on indefinite Finsler metrics, usually 
more implicitly than explicitly. The weakness of this approach is in the fact that there is no straightforward way of 
defining light rays in this setting. Asanov suggests a notion of light rays (see Chapter 7 in [2]) that depends on the 
choice of an auxiliary vector field. As the physical meaning of this auxiliary vector field is obscure, we do not think 
that this definition of light rays is satisfactory, from a physical point of view. Therefore, as we want to consider the 
equation of motion of light rays, we find Asanov's definition of Finsler structures inappropriate for the purpose of this 
paper. 

Fortunately, there is an alternative definition. Finsler metrics of Lorentzian signature were considered by Beem 
[4] in a way that is free from the above-mentioned drawbacks. In Beem's formalism there is no analogue of the 
Finsler function F; the Finsler structure is rather given directly in terms of the Lagrangian C{x, x), which should be 
sufficiently smooth (Beem requires it to be of class C 4 ) and real- valued for all x and all x ^ 0, it should be positively 
homogeneous of degree two, 

£(x, kx) = k 2 £(x, x) for all k > , (3) 

and the Finsler metric (JTJ) should be non-degenerate with Lorentzian signature for all x ^ 0. The non-degeneracy 
condition guarantees that the Euler-Lagrange equations © admit a unique solution to any initial condition (x(0), x(0)) 
with i(0) ^ 0. These solutions are the affinely parametrised Finsler geodesies which are well-defined for timelike 
(£(x,x) < 0), lightlike (£(x,x) = 0) and spacelike (£(x, x) > 0) tangent vectors x / 0. So in Beem's setting light 
rays can be unambiguously defined as lightlike geodesies, just as in standard general relativity. This is the reason 
why we consider, in this paper, a Finsler structure in the sense of Beem. Actually, for reasons that will become clear 
soon, we find it necessary to generalise Beem's definition a little bit: Our Lagrangian £ will not be smooth (and not 
even C ) at all (x,x) with x ^ 0; the second derivative of £ will give undetermined expressions on a set of measure 
zero. However, there will still be a unique solution curve (geodesic) through each point (x, x) with i ^ 0. 
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Another interesting generalisation of Beem's definition was brought forward recently by Pfeifer and Wohlfarth 
[171 118] . The main idea of their work is to allow for Lagrangians C that are homogeneous of any degree. However, as 
they still assume C to be "smooth" at all (x,x) with x =fi 0, their generalisation is of no advantage for our purpose, 
although it might be fruitful for other applications. 

Guided by Beem's definition [4] we define a Finsler spacetime in the following way. 

Definition 1. A Finsler spacetime is a 4-dimensional manifold M with a Lagrangian function C that satisfies the 
following properties: 

(a) £ is a real-valued function on the tangent bundle TM minus the zero section, i.e., C(x, x) is defined for all (x, x) 
with x =/= 0. 

(b) C is positively homogeneous of degree two with respect to x, i.e., eq. ([3]) holds. 

(c) The Finsler metric ([T]) is well-defined and has Lorentzian signature ( — h ++) for almost all (x, x) with i ^ 0. 
(As usual, "almost all" means "up to a set of measure zero".) 

(d) The Euler-Lagrange equations @ admit a unique solution for every initial condition (x, x) with x ^ 0; at points 
where the Finsler metric is not well-defined this solution is to be constructed by continuous extension. 

On a Finsler spacetime, we represent points in M by their coordinates x — (x° , x , x 2 , x 3 ) and points in the fibre 
T X M of the tangent bundle by their induced coordinates x = {x , x , i 2 , i 3 ). We use Einstein's summation convention 
for greek indices taking values 0,1,2,3 and for latin indices taking values 1,2,3. 

Note that the homogeneity condition (b) of the Lagrangian implies that the Finsler metric is positively homoge- 
neous of degree zero, 

g^vix, kx) = 9nv{x, x) for all k > , (4) 
and that the Lagrangian can be written in terms of the Finsler metric as 

£(x,x) = \^g liv {x 1 x)x^x v . (5) 

With the help of the Lagrangian we classify non-zero tangent vectors as timelike (C(x, x) < 0), lightlike (C(x, x) = 0) 
or spacelike (C(x, x) > 0). We call the solutions to the Euler-Lagrange equations @ the affinely parametrised Finsler 
geodesies. Again by the homogeneity condition (b) of the Lagrangian, C{x, x) is a constant of motion; hence Finsler 
geodesies can be classified as timelike, lightlike or spacelike. We interpret the timelike geodesies as freely falling 
particles and the lightlike geodesies as light rays. This interpretation is in agreement with the idea that the (Finsler) 
spacetime geometry tells freely falling particles and light rays how to move, i.e., that no additional mathematical 
structures enter into the equations of motion for freely falling particles and light rays. We have already mentioned 
that some authors disagree with this hypothesis, as far as light rays are concerned. As the interpretation of lightlikc 
Finsler geodesies as light rays is crucial for our work, some additional justification is given in the Appendix. 

In this paper we want to consider a special class of Finsler spacetimes that will serve us as a model for the 
gravitational field around the Sun. We shall assume that this gravitational field is static and spherically symmetric. 
In order to make these notions precise we have to recall that symmetries of Finsler metrics are described in terms of 
(Finsler generalisations of) Killing vector fields. By definition, a vector field V = V^d/dx 11 on a Finsler spacetime 
M is a Killing vector field if and only if its flow, if lifted to TM, leaves the Lagrangian C invariant. This condition 
can be rewritten in terms of the Finsler metric as 

8g pa dV* dg pa dV T dv T 

Here the depend on x only, whereas the g^ v depend on x and x. The Finslerian Killing equation ([B]) is known 
since the early days of Finsler geometry, see Knebelman [12] . 
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In the standard formalism of general relativity, one defines a spacetime as stationary if it admits a timelike Killing 
vector field V and as static if, in addition, this timelike Killing vector field V is orthogonal to hypersurfaces. If we 
exclude global pathologies (such as, e.g., the case that the quotient space M/V fails to be a Hausdorff manifold), the 
latter condition implies that the spacetime is a warped product of a 3-dimensional manifold with a (positive definite) 
Riemannian metric and the real line with a negative definite metric. We can use this property as the definition of 
staticity for Finsler spacetimes. 

Definition 2. A Finsler spacetime (M,C) is static if M is diffeomorphic to a product, M ~ R x N, and C is of the 
form 

C(x x ,x 2 ,x 3 ^i^ 1 ,x 2 ,x 3 ) = - (^gttix 1 , x 2 , x 3 ) i 2 + g ij (x 1 ,x 2 ,x 3 ,x 1 ,x 2 ,x 3 )x i x :i ^ , (7) 

where t runs over R and (a; 1 , a; 2 , x 3 ) are coordinates on N; the temporal metric coefficient g ft (x 1 , x 2 , x 3 ) must be 
negative and the spatial metric g%j{x , x , x , x , x , x ) must be positive definite. 

If the gij are independent of the x l , Definition [5] reduces to the definition of a static spacetime in the sense of 
general relativity. In any other case the limit of the for (i 1 , ± 2 , i 3 ) — s> (0,0,0) depends on the direction in which 
this limit is performed; this follows immediately from eq. Q. As a consequence, the Finsler metric fails to be 
well-defined on vectors tangent to the i-lines. This is the reason why, in part (c) of Definition [TJ the restriction to 
"almost all" non-zero tangent vectors was necessary to include proper Finsler Lagrangians of the form of eq. ([7]). 

We now add the condition of spherical symmetry. By definition, a Finsler spacetime is spherically symmetric if 
it admits a 3-dimensional algebra of Killing vector fields that generate the rotation group SO(3) such that each of 
its orbits is diffeomorphic to the 2-sphere S 2 . For a static Finsler spacetime as given in eq. ([7]). spherical symmetry 
means that we can choose the spatial coordinates as x 1 = r, x 2 = x 3 = ip, where r labels the group orbits and ~d 
and (p are standard coordinates on S 2 , and that then g t t depends on r only and the spatial part gijX l x 3 depends on 
r, r and ij 2 + sin 2 ?? <p 2 only. For a derivation of the latter fact see McCarthy and Rutz [TBI [25] . 



3 A class of spherically symmetric and static Finsler spacetimes 

As the class of all spherically symmetric and static Finsler spacetimes is too big, we make a more special ansatz for 
our model of the Solar system. We assume that the Lagrangian £ is of the form 

2C= (ha + c 2 ip )i 2 + ({hijhki + ip i j k i)x i x j x k x l \ 2 . (8) 

Here 

h tt dt 2 + hijdx^x 3 = h tt dt 2 + h rr dr 2 + r 2 (sin 2 tf dip 2 + dd 2 ) (9) 

is a spherically symmetric and static Lorentzian metric. In this section and in the following one, h tt and h rr are 
arbitrary functions of r, but later they will be specified to be the Schwarzschild metric coefficients, 

c2 2 /i 2GM , 

h ^~K r =- C (10) 

where c is the speed of light, G is the gravitational constant, and M is the mass of the gravitating body. The spatial 
perturbation ipiju is spherically symmetric and independent of t, 

ip ijM x i x j x k x l = ^i(r)r 4 + ^2(r)rV 2 (sm 2 i?^ 2 + i3 2 ) +^ 3 (r)r 4 (sm 2 i?^ 2 + '!? 2 ) 2 (If) 

and the time perturbation ipQ is a function of r only. 

Actually, ansatz ([8]) is less special than it might appear. The fourth-order term ^ijkii l x 3 x k x l can be viewed as 
the leading order term in a general Finsler power-law perturbation of the spatial part of the metric. (We do not want 
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to consider a third-order term because it would violate the symmetry under spatial inversions x l — > —x 1 .) For this 
reason, we consider the Lagrangian ([5]) as a natural choice for our purpose. 

In the following we refer to the dimensionless quantities t^A{f) as to the "perturbation functions", A — 0, 1,2,3. 
Throughout this paper, we assume that the perturbation functions depend differentiably on r and are so small that 
we may linearise all equations with respect to the iPa(t) and their derivatives tp' A (r). Differentiability and smallness 
of the guarantee that the Lagrangian C{x,x) is real- valued, and the Finsler metric ((TJ is non-degenerate with 

Lorentzian signature for almost all (x, x) with x ^ . The only points where this condition is violated are the points 
where the spatial velocity components are all zero, (i; 1 , a; 2 , a; 3 ) = (0,0,0), but t ^ 0. At these points, the Finsler 
metric gives undetermined expressions. We will see in the next section that, even through these points, the solutions 
to the Euler-Lagrange equations are uniquely determined by continuous extension, i.e., that our ansatz gives indeed 
a Finsler spacetime in the sense of Definition [T] 

We are still free to transform the radial coordinate. We can remove this freedom, thereby reducing the number of 
perturbation functions from four to three. In the unperturbed (Schwarzschild) spacetime, r is an "area coordinate" , 
i.e., the area of the sphere at r is given by Airr 2 . We can fix the radial coordinate by requiring that r has the same 
geometric meaning in the perturbed spacetime. From equations © and (|11[) we read that, in the perturbed spacetime, 
the sphere at r has area 47rr 2 (l -I- ^3). Hence, the desired condition is satisfied if we allow only perturbations with 
xp3 = 0. We are then left with three perturbation functions ipo, ipi and 02, and the Lagrangian §E§ reads 



2C = {h tt +c 2 ^)i 2 + (^h 2 r + ?Air 2 +r 2 (sin 2 tf <^ 2 -H9 2 )) 



\ ( y/h? r + ^if 2 + r 2 (sin 2 tf <p 2 + tf 2 )) ' 

From this expression we read that C is the Lagrangian of a pseudo-Riemannian metric if and only if 



(2h rr - 2^/h 2 r + ip 1 + *p 2 )r 2 f 2 (sin 2 tf ip 2 + d 2 ) 



(12) 



2h„ - 2^Jh 2 r + + ^2=0. (13) 



In this case the equations of motion can be investigated in terms of the standard PPN formalism. If (113)) does not 
hold, we have a proper Finsler geometry and the PPN formalism does not apply. We might say that the left-hand 
side of eq. (|13|) measures the "Finslerity" of our perturbed spacetime. 



4 Equations of motion 

We now discuss the solutions to the Euler-Lagrange equations © in our class of spherically symmetric and static 
Finsler spacetimes, i.e., the affinely parametrised Finsler geodesies. We restrict to timelike (C < 0) and lightlikc 
(C = 0) geodesies, which are to be interpreted as freely falling particles and as light rays, respectively. For timelike 
geodesies we can fix the parametrisation by requiring 2C = — c 2 ; then the affine parameter is equal to Finsler proper 
time r. 

By symmetry, it suffices to consider particles and light rays in the equatorial plane 1!) = tt/2. Then the linearised 
version of the Lagrangian (TT21 reads 

2C = (1 + Mh tt i 2 + (1 + ^)h rr r 2 + r V + t K f^f. 2 • (14) 

Here we have introduced, for notational convenience, modified perturbation functions 

c 2 ip ipt ip 2 h rr - -01 



h tt ' r 2h 2 r ' ^ z 2h 2 r 



(15) 



In terms of these modified perturbation functions, and after linearisation, the "non-Finsler condition" (fT5]) simply 
reads (f>2 — 0. Hence, in the linearised setting the "Finslerity" of our perturbed spacetime is measured just by fo. 
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By equation (|14[) . each of the perturbation functions tpQ, <pi and cf>2 has an obvious interpretation: <f>Q perturbs the 
time function t, <f>\ perturbs the radial length measurement and <p2 introduces a spatial anisotropy which is a genuine 
Finsler feature. Circular motion (f = 0) feels only </>o while radial motion (tp = 0) feels <po and <p±; the "Finslerity" 
<j>2 is felt only by motion that is neither circular nor radial. 

Equation (|14[) is the form of the Lagrangian on which all our following results are based. We will now derive the 
equations of motion. 

In addition to the constant of motion 



£ = 



for freely falling particles 



£ = for light, 

the t and tp components of the Euler-Lagrange equations give two more constants of motion E and L, 

dC 

— E = —r = (1 + 4>o) h tt t , 
at 



(16) 
(17) 

(18) 



L 



dC 

dtp 



r 2 tp 



1 



r 2 yj2j 



(19) 



The three constants of motion £, E and L give us three equations that determine the geodesies. From these three 
equations we read that, by continuity, there is a unique geodesic even for initial conditions f(0) = 0, tp(0) = and 
t(0) ^ 0, for which the Euler-Lagrange equations yield undetermined expressions, namely a curve with tp = constant. 
This completes the proof that our Lagrangian defines a Finsler spacetime in the sense of Definition [1] 

To within our linear approximation, the three conservation equations (fl4l . (|18l) and (IT91 can be solved for t, tp 
and r 2 , 



2£ 



<P 



htt 



t 



-E 
htt 



1 



L 



2£ 



htt r 



( E l 
\ htt 



J 



2C-*- 
h u 



2E< 



V htt 



(b E 2 
htt.hr, 



(20) 



(21) 



(22) 



From these three equations we find 



Equations (|23[) and (1241) determine the trajectories if parametrised by coordinate time t. If we are interested only in 
the geometrical shape of the trajectory, but not in its parametrisation, we may use the equation 



dip 



2C 



ii>rr 



AC- 



2E 2 3L 2 



\ 



V 



E 2 
htt 



2^4 



0o E 2 r' 

hrrhft.L 2 



(25) 



5 Circular orbits 



For a particle (2£ = — c 2 ) on a circular orbit, the equations dr/dip — and d 2 r/dtp 2 — must hold. By equation 
(|2"5")) . these two conditions are equivalent to 



htt 



I? 



boE 2 



0, 



hti U 



2L 2 



E- 



\htt 



With E 2 and L 2 determined this way, equation ((23]) yields 



/dip_Y 



Kt 

2r 



{ 



0o h' u 

h 2 
'hi 



t>oh t 



} 



After inserting the Schwarzschild metric (ITU)) , we find 



sdA£\ 2 _ GM 
Vdt) 



2GM 



(26) 
(27) 

(28) 
(29) 



If we denote the period by T, we have df/dt = 2tt/T, and fl2"9")) gives a generalisation of the third Kepler law for 
circular orbits, 



r 



{ 



1 - 



»o 



1 - 



2GM\ 



2GM 

In the unperturbed Schwarzschild spacetime, the Kepler law 

r 3 _ GM 
T 2 ~ 4tt 2 



))'} 



GM 
4vr 2 



(30) 



(31) 



coincides with the Newtonian Kepler law, as is well known. 

From an experimentalist's point of view, one may use the unperturbed Kepler law pip as an operational definition 
of GM. According to standard general relativity this would lead to a constant value GM. With our perturbation, it 
would lead to an r-dependent value GM{r) that is related to the constant GM value from general relativity by 



GM(r) 



GM\ 1 



2GM 



K-^))'} 



(32) 



We will now discuss the bounds imposed on (f>o(r) by (|32)) . 

From observations, the gravitational constant G is known, at present, only up to a relative uncertainty of approx- 
imatively 10 -4 . However, our knowledge of the product GM, where M denotes the Solar mass, is much better. The 
most recent value, taken from the webpage of the Jet Propulsion Laboratory http : //ssd. jpl .nasa.gov/?constants, 
is 

GM = 1.32712440018 x 10 20 mV 2 ± 8 x 10 9 mV 2 . (33) 
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More specifically, we get a value for GM = 4n 2 a 3 /T 2 from the observed values of the semi- major axis a and of the 
period T for each individual planet. As the periods are known with a higher accuracy than the semi-major axes we 

can write 

GM(a) - GM 



GM 



< 3e 



(34) 



where e is the accuracy with which the semi-major axis a is known. The values of e for the eight planets are shown 
in Tabled! 





Mercury 


Venus 


Earth 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


a/(10 10 m) 


5.79 


10.8 


15.0 


22.8 


77.9 


143 


288 


450 


Aa/m 


0.11 


0.33 


0.15 


0.66 


640 


4200 


3.8 xlO 4 


4.8 xlO 5 


e 


1.9X1CT 11 


3.1 xlO" 11 


l.OxlO- 11 


2.9X10" 11 


8.2xl0~ 9 


2.9xl0~ 8 


1.3xl0~ 7 


1.1 xl0~ 6 



Table 1: The first row of this table gives the semi-major axis a of each planet and the second row gives the formal 
standard deviation Aa of a. We have taken the values for Aa from Table 4 in Pitjeva [TH] who reported on data 
determined by the Institute of Applied Astronomy of the Russian Academy of Sciences from observations made 
between 1913 and 2003. In the third row we calculated the accuracy e of the semi-major axis for each planet, 
assuming that the real error may be one order of magnitude bigger than the formal standard deviation, e = 10Aa/a. 



In Section [5] below we will consider non-circular orbits. We will see that then </>i and 4> 2 do have an effect. For the 
time being we will be satisfied with a rough order-of-magnitudc estimate given by replacing each of the planets with 
a hypothetical planet that moves on a circular orbit with r = a. Then we can compare Q32p with (|34[) and conclude 
that 



(*>( 



c r / 



2 GM 
~ — tt^t- 3e 



Integration from t\ to r 2 yields 

fo(ra)(l 



2 GM 
c 2 r 2 



)-Mn)( 



l 



2GM\ 
c 2 ri J 



< 



2 GM | r 2 — n| 



3e r 



(35) 



(36) 



where e max is the maximal uncertainty between r\ and r-2- As 2GM/(c 2 r) varies from 10 8 near the Mercury orbit 
to 10 -10 near the Neptune orbit, (|36|) implies 



ri 



0o(r 2 ) - fain) 



1*2 



< 10 



-16 



(37) 



for all radii r\ and r 2 between the Mercury orbit and the Neptune orbit. If both r% and r 2 are between the Mercury 
orbit and the Mars orbit, the bound is even three orders of magnitude smaller. In particular, 



r \Mr)\ £ 10' 



- 16 



(38) 



between Mercury and Neptune, and even three orders of magnitude smaller between Mercury and Mars. Note that 
no assumption on monotonicity of the function <f>g was needed for this result. 

In analogy to the PPN formalism, one could assume that the perturbation functions are of the form 



2GM 



Al 



o 



A = 0,1, 2 



(39) 
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with constants 4>ai- We have included the Schwarzschild radius 2GM/c 2 in this ansatz to make the 4>ai dimensionless. 
Here we mean by GM the constant of Nature given by the fixed numerical value after the equality sign in just 
as we mean by c the constant of Nature given by the numerical value 299 792 458 m/s. If we accept the ansatz (|59l 
and if we assume that, as a reasonable first approximation, the terms of second and higher order can be neglected, 
the inequality (|35l) yields 

W(i-^)s».- w 

With the best value of e taken from Table [T] we find 

1 0oi I £ 3x HT 11 . (41) 



6 Radial free fall 

In this section we consider a freely falling particle (2C — — c 2 ) that moves in the radial direction, i.e. (p = 0. By (TK))) , 
the latter condition is equivalent to L = 0, so (|24p simplifies to 



2 p2 s A-^h 2 


\dtJ h rr E 2 \ ntt 



fdr\ 2 ft| / 2 £ 2 \, s 0oc 2 ^tt /.ox 



This expression gives the particle's velocity as it is measured by static observers with clocks that show coordinate 
time t. If we want to consider the same observers with clocks that show (Finsler) proper time r, we have to use the 
relation 

c 2 dr 2 = - {l + cj) Q )httdt 2 . (43) 

Then (|4"2"|) can be rewritten as 

/dr\2 c 2 h ft r/ , E 2 



JL{( c * + ^) (Wl)+c ^}. (44) 



KdrJ E 2 h r 

As a first application of this equation we want to discuss the acceleration of a particle from rest. From (1441) we read 
that, at a point where dr/dr = 0, the equation 

E 2 

c 2 + — + c 2 0o = (45) 
htt 

must hold. By differentiating (|44|) with respect to r, and inserting (1451) afterwards, we find the following expression 
for the acceleration from rest. 

^ d 2 r <^Kt | ^ + 0' ^" } (46) 

dr 2 huh rr I h' tt ) 

With the Schwarzschild metric (ITU1) this can be rewritten as 

d 2 r GM r 1 ,. c 2 r 2 ~i 

^2 — + (47) 



In the unperturbed Schwarzschild spacetime (I47D reduces to 

d 2 r GM 
dr 2 r 2 



(48) 



As an alternative to the method discussed in Section [SJ an experimentalist could use (|48| as an operational definition 
of GM. According to standard general relativity, this would lead to the same constant value GM as the method of 
Section [S] With our perturbation, however, it would lead to an r-dependent value 

2 2 

GM{r) = GM{ 1 - Mr) ~ <P'o(r) r + Mr)^ } (49) 
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which is different from the GM(r) of (I3"2"j) . 



GM(r) - GM(r) = GM { <f> (r) + 0i(r) } . 



(50) 



We see that the perturbation functions 4»q and can be determined, in principle, by observing circular orbits and 
radial acceleration from rest. With the bounds on cf>o we have found from the observation of circular orbits in Section 
[5J we can now discuss the bounds on </>i that result from measurements of free-fall accelerations. Again, we are 
satisfied with a rough order-of-magnitude estimate. We can then say that, with the Pioneer anomaly explained 
as a thermal recoil effect, all observations of radial accelerations up to the Neptune orbit are in agreement with 
General Relativity to within the order of the anomalous Pioneer acceleration of 9 x 10~ 10 m/s 2 . At the Neptune orbit, 
GM/r 2 » 6 x 10 _3 m/s 2 . As a consequence, (|4TJf suggests that 



Mr) 



2 

c r 
GM 



r$(r) 1 



2 GM 



Using \a\ - |b| < ||a| - |fe|| <\a- b\, we find 

\Mr) 

From Section [5] we know that 



< 



■ r #,(»■) (l - 



2GM\ 



< 



9 x 10 



-10 



6 x 10" 3 



+ 1.5 x 10 



-7 



(51) 



(52) 



2 GM 



GM 

r 4>' Q (r) iS 10 -6 between the Mercury orbit and the Neptune orbit, hence 

\Mr)\ & 10~ 6 (53) 



for all r in this range. 

If we assume that the perturbation functions have a fall-off behaviour according to (|39[) . and if we neglect terms 
of second and higher order, (|52p implies 



2GM 



< 2 



901 



1 - 



2GM\ 



') 



1.5 x 10" 



Evaluating at the Mercury orbit, 



2GM 



5 x 10 8 , and using (|H|) yields 



< 3 



(54) 



(55) 



This is much less restrictive than the bound on 0qi we na d found before. Note, however, that a value of <fin of the 
order of unity gives only a small correction to h rr of the Schwarzschild metric, because our ansatz (|39p involves the 
Schwarzschild radius. 



7 Effects on the paths of light rays 

In this section we want to calculate the effect of our Finsler perturbation on the worldline of a light ray that comes 
in from a source at radial coordinate rg, passes the Sun at a minimal value r m of the radial coordinate, and goes out 
again to an observer at radial coordinate ro- 

To that end we have to evaluate (|24l) and (|2"5|) with the Schwarzschild metric coefficients (|TU)) and C — 0. For 
notational convenience, we will use the abbreviation 

, / 2GM \ , , 

PW=r- 9 (l-- w -J (56) 
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throughout. With this abbreviation, and C — 0, equation (|25l) can be rewritten in the following form. 

( dr \ 2 4/^ \U . , 9 . o, 2 ^\ 0O^ 4 ^ 2 



We hrst determine how the constants of motion E and L depend on the minimum radius r m . If we insert the value 
r = r m into the right-hand side of (|57[) . we must get zero. If we solve the resulting equation for E 2 / X 2 , we find 



L 2 



c 2 p(r m ){l + Mr m )} ■ (58) 



Inserting this value into ([57)) yields 



where 
and 



2 = A (r){l-a(r)) (59) 

A(r) = r A (p(r m ) - p(r)) (60) 

t\ Af\ A f \ 2 P( r ™) ~ 3 ^( r ) (At \ A t \\ P( r ") fdt\ 

a(r) = </>i (r) - M r ) 7 — s {M r m) ~ Mr)) —, — ; rr ■ (61) 

P(r m ) P(r m ) - p{r) 



An analogous calculation puts (pMf into the form 

. dr 
~dt 



& =B(r)(l-m) (62) 



where 
and 

0(r) = Mr) ~ Mr)p(r f {rm ] ~ ^ - Mr) - (Mr m ) - Mr)) , ^ 7 ; ■ (64) 
Note that the perturbations a(r) and (3(r) depend not only on <po( r ), 4>i( r ) an d ^(O but also on <po(r m ), because of 



7.1 Light deflection 
From (1591) we find 



a(r) \ dr ,„_. 

' 1 + 4 r)7m (65) 



and integration yields the deflection angle A<p, 



* + / rS +Dfl + ^)^- (66) 



2 V ^Mr) 

If we denote by A</?o the deflection angle in the unperturbed Schwarzschild spacctimc, 

/ r rs r° \ dr 

* + Am = U„ + X.W (67) 

the deflection angle in the perturbed spacetime reads 

( f rs f r ° \ <*(r) dr 
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Using the mean value theorem, this result can be rewritten as 



= A(^ + ( 7T + A^o ) , (69) 



where f is some radius value between r m and max(rs,ro)- 

We want to evaluate these equations for the case that the source is a distant star (rg — > oo), the observer is on 
the Earth (ro = 1 AU), and the light ray is grazing the surface of the Sun (r m = 0.0046 AU). Then (j57]) gives the 
well-known deflection angle of 

Aipo = 1-75" = 8.48 x 10" 6 rad (70) 
and present day observations (see Will [29] , Section 3.4) require 

1 / . _ . . . \ - n _4 



Aip - Atp = g ( - 1-7 ± 4.5 ) x 10~ 4 A^o • (71) 



Comparison of (|69| and (fTTj) gives a bound for the possible values of a(f), 



A95 - Atpo 



Av? 



A(^o 



7T + Al^Q 



|a(f)| = 2 — ^ ^— < 2 x 10- 9 . (72) 



If we assume that the perturbation functions <^a(?") have a fall-off behaviour according to eq. p9[) . and if we neglect 
terms of second and higher order, the integral in (|68p can be calculated numerically. For r$ — > 00, ro = 1 AU and 
r m = 0.0046 AU we find 

A^ - A^o = 4.2 x 1(TVh + 5.1 x 10 _11 ^ 2 i - 4.2 x 10~%i • (73) 

We see that fai contributes with a much smaller factor than the other two perturbations. This has its reason in the 
fact that in (|6ip the Finslerity (j>2 (r) comes with a factor that changes sign between r — r m and r — ro, therefore 
positive and negative contributions to the integral partly cancel out. Hence, light deflection is rather insensitive to 
the Finslerity of our spacetime model. Combining (|73[) with (|7ip , and using (170)) and (|41l) , gives the quite insignificant 
bound 

I 8.2 x lO 4 0n + </>2i I ^ 52. (74) 

7.2 Time delay of light rays 
From (|62p we find 

1 rlr 

^=(l + -W)-^=. (75) 

Integration of this equation yields the travel time. The difference to the Newtonian travel time tjy is, by definition, 
the time delay St, 



t N + St = 



(r + r) (i+ ^ w »7fc- <76) 



In the unperturbed Schwarzschild spacetime, the time delay Stg is given by 



, prs pro . J r 



hence 

^«o + fr + rm. (78) 



/3(r) di 
2"7W 
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Using the mean value theorem, this can be rewritten as 



St 



Sto + ^P(r)(t N + St ) 



(79) 



where f is some radius value between r m and max(rs,ro). 

Time delays have been measured with radar signals since the 1960s. In the beginning Mars, Mercury and Venus 
were used as passive reflectors. In this case the round-trip travel time for a signal from the Earth to the planet and 
back is two times the one-way travel time (fT6")l plus a correction taking the orbital motion of the Earth into account. 
Later time delay experiments used spacecraft. The most accurate experiment of this kind was made with radio signals 
sent to the Cassini spacecraft, with the result (see Will [25], Sec. 3.4) that 



St - Stp 
St Q 



= (2.1 ±2.3) x 10" 



(80) 



The measurement was made when Cassini was at a distance of 8.43 AU from the Sun, and the distance of closest 
approach r m was 1.6 Solar radii (= 0.0074 AU). This corresponds to St ~ 273 /zs and ijv ~ 4700 s. Hence we find, 
with (1791) and (|80[) . a very small bound for a certain linear combination of the perturbation functions at some radius 
value r between 0.0074 AU and 8.43 AU, 



|0(f)| 



St - St 




St 


St 




t N + St 



< 5.2 x 10 



-12 



(81) 



If only the leading-order terms in (|39p are taken into account, numerical calculation of the integral in (|78p yields (for 
r m = 0.0074 AU, r s = 8.43 AU and r Q = 1 AU) 



St - St = ( 6.6 x 10~ 5 </>n + 3.3 x 10~ 6 fai - 7.5 x 10" 5 <foi ) s . 
The left-hand side can be estimated with the help of (1501) . With Sto ~ 273 (is and (|4"Tj) we find 



(82) 



20 < 



< 3.6 x 10" 



(83) 



8 Effects on bound orbits 

In Section [5] we have seen that circular orbits are affected only by the coefficient <po , but not by (f>i and fa ■ In 
this section we consider non-circular bound orbits, and we will investigate how the perturbation functions influence 
Kepler's third law and the perihelion precession. 

We consider a massive particle (2C = —c 2 ) on a bound orbit, with minimum radius r± (perihelion) and maximum 
radius T2 (aphelion). We need to calculate how the constants of motion E and L depend on n and r^- To that end, 
we rewrite (|25[) for the Schwarzschild metric coefficients, with 2C = —c 2 and using again the abbreviation (|56[) . in the 
following form. 

d 2 r^(E 2 -c^p-L 2 c 2 p) , 3hL 2 c 2 p x ^gy 

For r = r\ and r = r2 the right-hand side of (IM1) has to vanish, 

E 2 - c A r 2 p( ri ) - L 2 c 2 p( ri ) + Mn)E 2 - , E 2 - c 4 r 2 p(r 2 ) - L 2 c 2 p(r 2 ) + M^E 2 - . (85) 
Solving for E 2 and L 2 yields 

e 2 = E 2 + e 2 ( MriMn) - MnWz) \ 
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- i _ T -i , 3. ( Mr2)-Mn ) 

c 2 \ p(n) -p{r 2 ) 



where 



c 4 p(riMr 2 )(r|-r?) c 2 (p(r 2 )r| - p{r x )r\) 



E o = — ^^-^-^^——Ll an d Lq 



p(ri)-p(r 2 ) ^ p(ri)-p(r 2 ) 

are the values of the constants of motion in the unperturbed Schwarzschild spacetime. Substitution of (|86)l and (|87| 
into (|84|). and using the identities E$ = L\(?p(r\) + c A p(ri)r\ — Llc 2 p(r 2 ) + c 4 p(r 2 )r%, gives the orbit equation in 
terms of Eq and Lq, 



2 P 



c 2 E, 



where 



and 



r(ri + r 2 — r) 
2GM \ ( 2GM 



2 / 2UM \ f ZUM \ 

cr(r 1 +r 2 )\l--^-j\l-^-j /r 2 (0 o (r 2 ) - o (r)) n(0o(»l) - MO) 



^(m-r.j/l-^ri + i + i)} 



2GM / 1 1 1 

— (- + - + - 

&■ \r\ r 2 r 

For the limiting case of a circular orbit, r = r± = r 2 , we find 



r 2 — r r% — r 



1 



Analogously, substitution of (|5r?)) and (I57|) into results in 



where 



2Gtf r 2 -r r-n 1 5— 1 — - + - 

V c z r I \ c z V ri To 

r 2 (r 1 +r 2 ) ( 1 - ) ( I 



c 2 ri J \ c 2 r 2 



(89) 



^_ r±(E 2 -c±r 2 p~L 2 c 2 p) , 3 2 Lge 

o{(0o(r 2 ) -0oW)(p( r i) r i ~p(r)r 2 ) - (</>o (ri) - Mr)) {v{r 2 )rl -p(r)r 2 )} 
L 2 (p(n) - p(r 2 )) (e 2 - c*r 2 p - L 2 c 2 p) 
After substituting an d Lq from (|88|) we find 

2 = C(r)(l- 7 (r)) (90) 



C(r) = ^-^-^) f ! _ ^(1 + 1 + 1)) (91) 
ri r 2 \ c z V n r 2 r / / 

. . / 2GM\ 
3(p 2 (r)rir 2 [ 1 — 

7 (r) = fc(r) - 2Mr) + 7 — —\ ry (92) 

2GM (rf + r% + rir 2 — rr\ - rr 2 ) \ 

c 2 nr 2 (ri +r 2 - r) J 



7 (r) = Mr) + Mr) + — ^ otH <#>(r) + U») ■ (93) 



(^) 2 =U(r)(l-*(r)) (94) 



2GM\ 2 ( 2GM ( 1 1 1 



D{r) = ^ 2GM \ r 2GAf x (95) 
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and 



«(r) = 7 (r) - 2*,(r) 2GM(r2 _ ri) + ( 96 ) 



2</> 2 (r)(r 2 -r) 2 (r-r 1 ) 2 (^l - — (- + - + -)J 

r 2 [ n + r 2 - r - ^„ (r\ + r\ + rir 2 - r(r\ + r 2 )) ] 
\ c z r x r 2 J 



For the limiting case of a circular orbit, r = r± = r 2 , we find 

*(r) - 7(r) = -to(r) - ^ (l - ) &(r) . (97) 

Equations (|93| and (|97| can be used as valid approximations for orbits whose eccentricity is not too big, where r is 
any value between the perihelion and the aphelion. 

8.1 Perihelion precession 
From (jnnj) and JM} we find 

^=(l+h(r))J= (98) 



dt 



(l+l S(r) ) . (99) 



Integrating these two equations over the orbit from one perihelion transit to the next, 

2?r + A<& 



T = 2 



( 1 + r<'->) Jft ( 10 »> 

i 1+ \^)7m- (101) 



gives the anomalistic period T (in terms of coordinate time t) and the angular advance A$ of the perihelion during 
this period. We denote the corresponding quantities in the unperturbed Schwarzschild spacetime by an index 0, 

r T2 dr 

2tt + A$ = 2 / -_, (102) 
r 2 dr 

To =2/ -_. (103) 

The precession rate of the perihelion, in radians per time, is uj — A$/T. In our linearised setting uj deviates from 
cj = A<5> /T according to 

oj-oj q _ 1 f 2 7 (r)dr 1 f 2 5(r)dr 



uj A$ J ri T J n y/Dffi 

With the mean-value theorem the last equation can be rewritten as 



U "° - lif) ( 2n ' 0-^, (105) 



too 2 V A$ J 2 



1G 



where f and f are some radius values between n and r%. 

For Mercury, the precession rate ujq according to general relativity is well-known to be 43 arcseconds per century. 
This corresponds to a precession angle (in radians) per revolution of A$o = 0.502 x 10 -6 . Present day observations 
confirm that the general-relativistic value is true, with a possible relative error of 10 -3 , see Will [55]; Section 3.5. 
Hence, (|105[) implies 

|6.2 x 10 6 7(f) - 0.5 (5(f) | < 10~ 3 , (106) 

which is a restrictive bound on 7(f) for r on the Mercury orbit. 

If we take only leading-order terms of (j3"9")l into account, the integrals in (|104p can be numerically calculated. With 
the Mercury values (A$ = 0.502 x 10~ 6 , T = 87.969 d, n = 46 001200 km and r 2 = 69 816 900km) we find 

- = 3.3 x 10 _1 0ii + 3.1 x 10 _1 02i + 2.5 x 10~ 8 O i . (107) 

Wo 

The observational fact that the left-hand side is bounded by 10~ 3 , toghether with (|4"Tj) . implies that 

I 0ii + 9.4 x lO _1 02i I ^ 3.0 x 10~ 3 . (108) 
In combination with (|83[) this gives us a bound for 02 i> 

I 02i I 3.6 x 1CT 3 , (109) 

which means that the Finslerity is bounded by 

I 2 (r) I < 1.8 X 10~ 10 (110) 
everywhere in the Solar system beyond the Mercury orbit. 



9 Conclusions 

In this paper we have considered a class of spherically symmetric and static Finsler spacetimes which are small 
perturbations of the Schwarzschild metric. After fixing the ambiguity in the choice of the radial coordinate by 
requiring that a sphere at coordinate r has area 47rr 2 , the perturbed metric is characterised by three functions 0o(r), 
0i (r) and 02 (r) which we called the "perturbation functions" . It was our main goal to determine the bounds which 
are imposed on these perturbation functions by observations in the Solar system. In this way we have provided a 
framework for testing if a certain Finsler modification of general relativity is in agreement with experimental facts. 

We have been careful to set up the formalism in such a way that not only freely falling particles but also light 
rays are unambiguously defined as Finsler geodesies. We feel that this is a major advantage in comparison to several 
other Finsler approaches where the definition of light rays is questionable. Having both freely falling particles and 
light rays at our disposal is essential because these are the tools needed for the experiments discussed. 

The formalism presented here is meant as an analogue of the PPN formalism. From a methodological point of 
view, there are two differences. First, our formalism is post-Schwarzschild rather than post-Newtonian. This is, of 
course, motivated by the fact that we wanted to concentrate on the possible Finsler deviations from standard general 
relativity. Second, we chose for the radial coordinate the area coordinate, whereas in the standard PPN formalism 
one chooses the isotropic radial coordinate. This is a necessary deviation from the standard PPN formalism because 
the isotropic radial coordinate does not exist in a proper Finsler spacetime. 

Our approach is purely kinematical, i.e., no field equation is used. Hence, one can use it for testing the validity of 
solutions to any Finslerian field equation, provided that the solutions belong to the class considered in this paper. 

Here we have discussed only tests where all objects moving in the field of the Sun can be treated as test particles. 
One could set up a Finsler geometry model for more complicated situations, e.g. for the motion of the Earth in the 
combined gravitational field of the Sun and the Moon. This would make more sensitive tests possible, in particular 
using the very precise Lunar Laser Ranging measurements. This is planned to be done in future work. 
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Appendix 



For our work it was crucial that the timelike Finsler geodesies were interpreted as freely falling particles and the 
lightlike geodesies were interpreted as light rays. We feel that this is the most natural interpetation, if one believes 
in the possibility that our real world carries a Finsler spacetime structure. As far as the light rays are concerned, 
the geodesic hypothesis can be further justified by deriving the lightlike geodesies as the bicharacteristic curves of 
appropriately generalised Maxwell equations. In this appendix we will outline how this can be done. As a detailed 
treatment would require a separate paper, we will only sketch the line of thought, just enough to convince the reader 
that the geodesic hypothesis for light rays is, indeed, well motivated. 

In Definition[T]we have defined Finsler spacetimes in terms of a Lagrangian C(x, x). We can switch to a Hamiltonian 
formulation by introducing canonical momenta 

dC{x,x) 

Pm = gifl = 9fj,u{x, x)x (111) 

and the Hamiltonian 

H(x,p) — p^x^ 1 £{x 7 x) = -g p,u (x,p)p tl p u (112) 

where g^(x,p) is defined through 

g^(x,p)g ua (x,x)=5Z. (113) 

In (|112l) and (|113l) . i M must be expressed as a function of x and p with the help of (| 11 1[) . For later convenience, we 
write 

H^p) = «^ . (114) 

Note that, because jC(x,x) is assumed to be homogenoeus of degree two with respect to the i M , the Hamiltonian 
T-L(x,p) is homogeneous of degree two with respect to the p^, hence 

p^(x,p)=2U(x,p). (115) 

The lightlike Finsler geodesies are the solutions to Hamilton's equations with H(x,p) — 0. It is our goal to demonstrate 
that these are the bicharacteristic curves of appropriately generalised Maxwell equations. 

In the case of a pseudo-Riemannian metric, where the g^ depend on x only and not on p, the source-free vacuum 
Maxwell equations can be written as 

d»F w {x) + d v F^{x) + d.F^x) = (116) 

(n^(x,d) + ... )F fiU = 0. (117) 

Here F^ v is the electromagnetic field tensor and <9 M means partial derivative with respect to x^. In (j 1 1 T[) we have 
written only the principal part. H M (a;, d) = g llcr (x)d a is a first-order differential operator and, thus, homogeneous of 
degree one with respect to the d a . The omitted terms, indicated by ellipses, involve the Christoffel symbols and no 
derivatives, so they are in particular homogeneous of degree zero with respect to the <9 M . 

It is very natural to postulate that Maxwell's equations take the same form of (|116[) and (|117[) on a Finsler 
spacetime. Now "H M (a;,p) is no longer a polynomial, but still homogeneous of degree one with respect to the p^, 
so W 1 {x,d) is no longer a differential operator but still a well-defined pseudo-differential operator. (For a detailed 
exposition of pseudo-differential operators see, e.g., Hormander [H].) The terms indicated by ellipses are necessary 
in (|117l) to make this equation coordinate-independent. Their special form, however, will not be relevant for the 
following argument. We only have to assume that they are homogeneous of degree one with respect to the d^, as in 
the pseudo-Riemannian case. 

We now apply the operator H^{x, d) to ([TT6| . With the help of (fTT7|) and (fTT5)l we find 

U{x,d)F ua + ... = (118) 
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where again only the principal part (i.e., the terms of highest degree of homogeneity) has been written out. This 
demonstrates that the field tensor satisfies a (pseudo-differential) Finslerian wave equation. The principal part 
determines the characteristic equation (or eikonal equation) 

H(x,dS) = 0. (119) 

It gives the characteristic surfaces S = const along which solutions F IJiV {x) to Maxwell's equations might have dis- 
continuities of their first derivatives. The characteristic equation has the form of a Hamilton- Jacobi equation with 
the Finsler Hamiltonian Ji(x,p). The bicharacteristic curves (or rays) are the corresponding solutions to Hamilton's 
equations, i.e., the lightlike Finsler geodesies. 

We have thus derived the result that light rays are lightlike Finsler geodesies from the assumption that Maxwell's 
equations on a Finsler spacetime take the form of (|116|) and (|1 1T|) . A full treatment of the subject would, of course, 
require to specify the omitted terms in (|117p and to discuss their implications for physics. This could be the subject 
of another paper. 

Acknowledgments 

CL and VP gratefully acknowledge support from the Deutsche Forschungsgemeinschaft within the Research Training 
Group 1620 "Models of Gravity" . Moreover, VP was financially supported by the German-Israeli Foundation, Grant 
1078/2009, and by the Deutsche Forschungsgemeinschaft, Grant LA905/10-1, during the course of this work. CL 
wishes to thank the cluster of excellence QUEST for support and Jiirgen Miiller for discussions. 

References 

[1] A. K. Aringazin and G. S. Asanov. Finslerian post-Riemannian corrections to the equations of geodesies. Gen. 
Relativ. Gravit., 17:1153-1163, 1985. 

[2] G. S. Asanov. Finsler geometry, relativity and gauge theories. Reidel, Dordrecht, 1985. 

[3] G. S. Asanov. Finslerian extension of Schwarzschild metric. Fortschritte der Physik, 40:667-693, 1992. 

[4] J. K. Beem. Indefinite Finsler spaces and timelike spaces. Can. J. Math., 22:1035-1039, 1970. 

[5] G. Y. Bogoslovsky. On a special-relativistic theory of anisotropic space-time. Dokl. Akad. Nauk SSSR, 213:1055- 
1058, 1973. 

[6] Z. Chang and X. Li. Modified Newton's gravity in Finsler space as a possible alternative to dark matter hypothesis. 
Phys. Lett., B 668:453-456, 2009. 

[7] A. A. Coley. Clocks and gravity. Gen. Relativ. Gravit., 14:1107-1114, 1982. 

[8] F. Francisco, O. Bertolami, P. J. S. Gil, and J. Paramos. Modelling the reflective thermal contribution to the 
acceleration of the Pioneer spacecraft. Phys. Lett., B 711:337-346, 2012. 

[9] G. Gibbons, J. Gomis, and C. N. Pope. General very special relativity is Finsler geometry. Phys. Rev., D 
76:081701, 2007. 

[10] F. Girelli, S. Liberati, and L. Sindoni. Planck-scale modified dispersion relations and Finsler geometry. Phys. 
Rev., D 75:064015, 2007. 

[11] L. Hormander. The analysis of linear partial differential operators III, IV. Springer, Heidelberg, 1985. 
[12] M. S. Knebelman. Collineations and motions in generalized spaces. Amer. J. Math., 51:527-564, 1929. 



19 



[13] V. A. Kostelecky. Riemann-Finsler geometry and Lorentz-violating kinematics. Phys. Lett., B 701:137-143, 2011. 

[14] C. Lammerzahl, D. Lorek, and H. Dittus. Confronting Finsler space-time with experiment. Gen. Relativ. Gravit., 
41:1345-1353, 2009. 

[15] X. Li and Z. Chang. A possible scenario of the Pioneer anomaly in the framework of Finsler geometry. Phys. 
Lett, B 692:1-3, 2010. 

[16] P. J. McCarthy and S. F. Rutz. The general four-dimensional spherically symmetric Finsler space. Gen. Relativ. 
Gravit, 25:589-602, 1993. 

[17] C. Pfeifer and M. N. R. Wohlfarth. Causal structure and electrodynamics on Finsler spacetimes. Phys. Rev., D 
84:044039, 2011. 

[18] C. Pfeifer and M. N. R. Wohlfarth. Finsler geometric extension of Einstein gravity. Phys. Rev., D 85:064009, 
2012. 

[19] E. V. Pitjeva. High-precision ephemerides of planets - EPM and determination of some astronomical constants. 
Solar System Research, 39:176-186, 2005. 

[20] B. Rievcrs and C. Lammerzahl. High precision thermal modeling of complex systems with application to the 
ffyby and Pioneer anomaly. Ann. Phys. (Leipzig), 523:439-449, 2011. 

[21] I. W. Roxburgh. Post-Newtonian limit of Finsler space theories of gravity and solar system tests. Gen. Relativ. 
Gravit, 24:419 431, 1992. 

[22] I. W. Roxburgh and R. K. Tavakol. Non-Riemann geometrizable effects in the gravitational one-body problem. 
Gen. Relativ. Gravit, 10:307-312, 1979. 

[23] H. Rund. The differential geometry of Finsler spaces. Springer, Berlin, 1959. 

[24] H. Rund and J. Beare. Variational properties of direction- dependent metric fields. University of South Africa, 
Cape Town, 1972. 

[25] S. Rutz. Symmetry in Finsler spaces. In D. Bao, S.-S. Chern, and Z. Shen, editors, Finsler geometry, volume 196 
of Contemporary Mathematics, pages 289-299, Providence, Rhode Island, USA, 1996. American Mathematical 
Society. 

[26] S. Rutz. A Finsler generalisation of Einstein's vacuum held equations. Gen. Relativ. Gravit., 25:1139-1153, 1998. 

[27] L. I. Schiff. On experimental tests of the general theory of relativity. Amer. J. Phys., 28:340-343, 1960. 

[28] S. G. Turyshcv, V. T. Toth, G. Kinsella, S.-C. Lee, S. M. Lok, and J. Ellis. Support for the thermal origin of 
the Pioneer anomaly. Phys. Rev. Lett, 108:241101, 2012. 

[29] C. M. Will. The confrontation between general relativity and experiment. Living Rev. Relativity, 9(3), 2006. 
http: / /www.livingreviews.org/lrr-2006-3. 



20 



